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Abstract. The stresses in a closed lipid membrane described by the Helfrich 
hamiltonian, quadratic in the extrinsic curvature, are identified using Noether's 
theorem. Three equations describe the conservation of the stress tensor: the 
normal projection is identified as the shape equation describing equilibrium 
configurations; the tangential projections are consistency conditions on the 
stresses which capture the fluid character of such membranes. The corresponding 
torque tensor is also identified. The use of the stress tensor as a basis for 
perturbation theory is discussed. The conservation laws are cast in terms of 
the forces and torques on closed curves. As an application, the first integral of 
the shape equation for axially symmetric configurations is derived by examining 
the forces which are balanced along circles of constant latitude. 



PACS numbers: 87.16. Dg, 46.70.Hg 
1. Introduction 

The Helfrich hamiltonian, quadratic in the extrinsic curvature, provides a remarkably 
robust description of the mechanical properties of lipid membranes (see, for example, 
|[], ||, ||, [|, 0], and the comprehensive review 0). This hamiltonian associates 
an energy penalty with bending. It is invariant with respect to deformations 
along the membrane itself. Equilibrium configurations will therefore satisfy a single 
'shape' equation corresponding to the extrema of the hamiltonian with respect to 
a deformation directed along the normal^, This equation is purely geometrical. 
Underpinning any non-trivial geometry however, there will always be stresses within 
the membrane. In this paper, we will focus on the distribution of these stresses, as well 
as the corresponding torques, and how the two are reflected in the geometry. Apart 
from its conceptual value, it might appear that the formulation of the problem this 
way represents a step in the wrong direction; after all, the shape equation describes 
the physics completely, and to obtain it, there is no need to know what the stresses 
are, much less the torques. It is surprising then, that by stepping back a little from 
the geometry this way, a potentially powerful as well as economical approach to the 
physics is also provided. We should remark that the idea itself in this context is not 
new, going back to the early work by Evans jl(|. More recently, it has also been 
exploited for axisymmetric geometries in to study fluid membrane tethers. 

TTo construct the stress tensor we appeal to Noether's theorem: the euclidean 
invariance of the Helfrich hamiltonian implies conservation laws and corresponding 
conserved charges. (For an analogous treatment of elastic curves in space see In 
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particular, translational invariance implies the proportionality of the divergence of the 
membrane stress tensor to any applied external forces; if the membrane is closed and 
the enclosed volume is fixed, this source will be the hydrostatic pressure enforcing the 
constraint. If the membrane is acted on by a localized external force, as is the case of 
micro-manipulation techniques 12, D|, this will appear as an additional distributional 
source on the membrane. 

Implicit in Noether's theorem is the recipe for the construction of the relevant 
conserved quantities. This is important because the corresponding Newtonian 
construction, while simple in principle, is not in practice. (See, however, |jl4|| , where 
the stresses and the torques arc derived for axisymmctric configurations.) 

The relationship between the three equations describing the conservation of the 
stress tensor and the (single) shape equation is very simple: the projection onto 
the normal of the former is the shape equation. The projection process necessarily 
dismantles the divergence form of these equations: the divergence of the normal stress 
picks up a source proportional to the local tangential stress, the 'extrinsic curvature 
cubed' nonlinearity in the shape equation. The fluid character of a lipid membrane 
is captured mathematically in the reparametrization invariance of the hamiltonian 
describing it. The two tangential projections of the conservation law (or its 'Bianchi 
identities' in the language of gauge theory) reflect this invariance. As such, these latter 
two equations are satisfied even when the shape equation is not. This structure itself 
is model independent. There is thus little or no cost involved in considering a general 
hamiltonian depending on the extrinsic curvature, in particular, models higher order 
in curvature. 

It is also possible to cast the conservation of the stress tensor as a more intuitive 
global statement relating the total force on any closed loop lying on the surface to the 
action of external forces on the area enclosed by the loop. This global statement of 
conservation is particularly useful when the membrane possesses some level of spatial 
symmetry. The first integral of the shape equation for axially symmetric configurations 
of both spherical and toroidal topology is obtained as an immediate consequence when 
the loop coincides with a circle of constant latitude. This approach to the problem 
contrasts favorably with 'hamiltonian' approaches jl5], [l6|]. We also identify a global 
constraint on the tangential stress over the closed membrane, which reproduces the 
known scaling identity for the Helfrich hamiltonian . 

The paper is organized as follows. In Section 2 we describe how the hamiltonian 
responds to a deformation of the membrane surface. In Section 3 it is shown how 
Noether's theorem may be applied to identify the conservation of the stress tensor 
as the conservation law associated with translation invariance of the two-dimensional 
surface. In this section, we also show how the conservation of the stress tensor can 
be projected to cast the shape equation in terms of derivatives of the stress tensor 
components. The special case of a soap bubble is used in Section 4 to illustrate our 
approach. In Section 5, we extend our considerations to the Helfrich hamiltonian and 
more generally to hamiltonians that depend on any power of the extrinsic curvature. 
We derive explicit expressions for the stress tensor. In Section 6, we consider rotational 
invariance and the definition of torque. We introduce an intrinsic torque which, added 
to the couple due to the stress tensor, is conserved. In Section 7, the global form of 
the conservation laws is examined; a non-trivial global relationship which involves the 
stress tensor itself is identified. In Section 8 we comment briefly on how the stress 
tensor might be exploited to develop perturbation theory. In Section 9, we specialize 
to axially symmetric configurations. We show how to obtain a relationship between 
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the components of the stress tensor that encodes the shape equation. We conclude 
with some brief remarks in Section 10. 



2. Response of the hamiltonian to surface deformations 

Let us consider an embedded surface £ in three-dimensional euclidean space. For 
simplicity, we suppose that this surface is closed. The surface is described locally by 
three functions x = X(£ a ), where x = (a; 1 , x 2 , a; 3 ) are for now cartesian coordinates 
on R 3 , £ a ,a = 1,2 may be any local coordinates on the surface. We introduce two 
tangent vectors e a which are defined by e Q = <9X/9£ a . The surface geometry is 
described completely by the induced metric j ab , and the extrinsic curvature K ab . The 
former is given by 

lab — e Q • e b . (1) 

We denote by V a the covariant derivative on the surface which is compatible with j ab . 
The extrinsic curvature K ab is given by 

K ab = -n • d a e b , (2) 

where n is the unit vector normal to the surface. We will indicate the trace of K ab by 
K = j ab K ab . The intrinsic and extrinsic geometries are related by the Gauss-Codazzi- 
Mainardi equations, 

K ab K ab -K 2 -U = Q, (3) 
V a K bc - V b K ac = , (4) 

where 1Z denotes the surface scalar curvature. 

Though our principal interest in is the case of a lipid membrane, we will consider 
any local hamiltonian, depending on the functions X, which is invariant under surface 
reparametrization; no energy penalty is associated with shearing deformations of the 
membrane. The surface hamiltonian is given by 

JT E [X] = / dA h ( 7o6 , K ab , V a K bc , ■■■). (5) 

where the scalar h is constructed locally from the geometry of the surface, and the 
infinitesimal area element is given by 

dA = d 2 ^, (6) 

with 7 = det7 {,. In particular, the Hclfrich hamiltonian is proportional to 

dA K 2 . (7) 

E 

If there is a spontaneous curvature Co, K is replaced by K — Co in Eq.(0). If the 
membrane is closed there may be a contribution to the total hamiltonian proportional 
to the enclosed volume V due to a constant pressure excess P there. In addition, there 
may be global constraints imposed on the membrane geometry. For example, one or 
more of the area, the integrated mean curvature, 



M= / dAK, (8) 



s 



or the enclosed volume V may be fixed. In Hclfrich's original model, both A and V 
are fixed. In Svetina and Zek's bilayer couple model, in addition to these M gets 
fixed |n| . These constraints are imposed by lagrange multipliers. We note that in the 
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refinement of these models, known as the area difference model, the hamiltonian is 
not expressible in the simple form (JsJ) : a non-local term proportional to (M — M ) 2 is 
added to the hamiltonian, where M is some constant fjj, H |(], In this model, 
deviations in the integrated mean curvature from Mq are penalized energetically. We 
will not consider this model further. 

To proceed, we need to know how this hamiltonian responds to an infinitesimal 
deformation of the embedding functions for the surface, X(£) — > X(£) + 5X(£). We 
decompose <5X into its parts tangential and normal to the surface, 

SX = $ a e Q + $n . (9) 

As an intermediate step, we evaluate how the quantities "/ ab and K ab respond to this 
deformation. The tricky bit is to do this in a way which does not depend on the 
particular reparametrization of the surface we choose. For the induced metric we have 

^lab = e Q • V b <5X + e b ■ V SX 

= 2 K ab $ + V a $ b + V b $ a . (10) 

To obtain the second line we have used the decomposition (Q), together with the 
Gauss- Weingarten equations, 

V a e b = - K ab n, (11) 
V a n = K ab e b . (12) 

Similarly, for the extrinsic curvature defined by (|2j), we find 

SK ab = - (<5n) • V Q V 6 X - n • V a V b 5X 

= - V Q V b $ + K ac K c b <S> + <S> c V c K ab + K ac V b <S> c + K bc V a <$> c . (13) 

Here the first term in the first line vanishes, because of the unit vector fact: 
(<5n) • V a V b X = —K ab (<5n) • n = 0. To obtain the third term in the second line 
we have used the Codazzi-Mainardi equation As expected, in both Eqs. (|i~0|), 
(^3|), the tangential deformation is the Lie derivative along the surface vector $ a . 

Once we know how the geometry changes under a surface deformation, we are 
in a position to vary any hamiltonian of the form (^|). However, we can exploit 
reparametrization invariance to simplify matters. The infinitesimal change in the 
hamiltonian can always be decomposed into its tangential and normal parts, 

5H = S\\H + 5 ± H . (14) 

Away from boundaries, the tangential deformation can be identified with a 
reparametrization of S since 8\\H is a boundary term. It is simple to show why 
this is so: we have that 8\\f = & a d a f for any scalar function /(£) defined on S; in 
addition, under a tangential deformation, the induced metric on S tranforms as a Lie 
derivative, as given by the tangential part of (|l^). Thus we have that 

S\\Vl = UVl ® a ) ■ (15) 

The identification of a conservation law will require us to isolate a divergence in 
the variation of the hamiltonian. Thus let us consider the behavior of the hamiltonian 
iJj] when it is restricted to some connected domain So with a boundary C. In the 
case of a surface with spherical topology, we will suppose without loss of generality 
that Eo is also simply connected so that C is a contractable closed curve. For toroidal 
or higher genus surfaces, one might also be interested in regions which are not simply 
connected with a non-contractable, disconnected boundary. For notational simplicity, 
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we will proceed as though C is also connected. For an arbitrary deformation of E we 
have 

SH So = [ d 2 H {(S^f) h + ^j (Sh)} . (16) 
A tangential deformation of the surface thus always results in a pure divergence, 



*||flEo= / d 2 tid a (^h1> a ) = [dshl a $ a , 



(17) 



where we have used Stokes theorem in the second equality, s is arclength along the 
boundary curve C, induced by its embedding in S , say £ a = Y a (s), and l a is the unit 
normal on C pointing out of So. 

A reparametrization of So can only move its boundary. Thus, if there is no 
boundary, 5\\H-z = 0. If we are interested in a stiff membrane, so that no local 
change in the area is permitted, then <!> a itself must respond to $ to maintain the 
constraint, = 0. Then, as follows from Eq. jlC|), V a $ a + K<& = 0. Clearly 

this implies no constraint for a euclidean motion. For general deformations, we only 
require that the tangential deformation results in a divergence. Thus none of our 
conclusions are modified. The enclosed volume also clearly does not change under a 
surface reparametrization. 

Whereas the tangential variation of the hamiltonian is simple, the normal 
variation is, in general, non-trivial. However, the latter can always be cast in the 
form 

5xHx = [ dA {8(h) $ + V a S a {<S>] } , (18) 

i.e. as a bulk part plus a pure divergence. Here 8(h) is the Euler-Lagrange derivative 
of h with respect to surface deformations, projected onto the normal n to the surface; 
S a is a linear differential operator on the surface which operates on the normal 
deformation $ as follows: 

S a [t>] = 5 ( Q 0) $ + S$V 6 $ + • • • . (19) 

To construct S' a [ < i>], integration by parts is used to collect in a pure divergence all 
surface gradients and higher derivatives of <I>. Such terms will show up in the normal 
deformation of K a \, and its derivatives. 

Summing the two independent variations we have 

6H So = f dA [8(h) $ + V Q Q a ] , (20) 

where the Noether charge is 

Q a = S a [$} + h <f> a . (21) 

Note that Q a is not unique; clearly Q a — > Q a + e afc Vb/, for some scalar density /, 
with e ab the surface Levi-Civita density, will leave Eq. ( po| ) unchanged. 

The variational principle restricted to normal deformations determines the 
equilibria. Suppose that the enclosed volume is fixed at some value Vq. The total 
hamiltonian 

H = H^-P{V- Vq) (22) 

is then stationary with respect to normal deformations of S when the Euler-Lagrange 
equation 

8(h) = P (23) 
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is satisfied. This is the 'shape' equation determining the equilibria of the membrane. 
We have used the fact that the interior volume varies as 



5±_V = dA $. (24) 
3. Stresses in the membrane 

To this point, we have considered arbitrary deformations of the functions X. We 
now ask what occurs if we subject the membrane to a translation or a rotation. Let 
us then consider an infinitesimal euclidean motion, = a + b x X, where a is an 
infinitesimal constant translation, and b represents an infinitesimal rotation. Let us 
focus on translations. Rotations will be considered in Section ^|. 

If we decompose 5X according to Eq.(^J) then for an infinitesimal spatial 
translation, <5X = a, we have $ = a n, and $ a = ae a . Substituting into Eq. (|20|) , the 
variation of the hamiltonian associated with this translation can be cast in the form, 



6Hv =a- / dA [£(h)n-V a f a ] 



(25) 



The surface vector f a is given by 

f a = -S a [n] - he a , (26) 
where S a [n] is defined by 

S a [<S>] =a-S a [n]. (27) 

The quantity f a describes the non-vanishing components of the stress tensor as 
appropriate for a two-dimensional system: it is both a spatial vector and a surface 
vector. The ambiguity in the Noether charge Q a translates into an ambiguity in the 
stress f a — > f a + 6 ab V{,W, for some arbitrary vector density W. 

While tangential deformations do not participate in the variational derivation of 
the shape equation, we see that they do contribute in a simple but essential way to 
the construction of the stress tensor. 

Suppose that Hy, is invariant under translations, 5Hy, = 0. Because So is 
arbitrary in Eq. (p5f) , the integrand vanishes pointwise so that 

V a f a = £(h)n. (28) 

If, in addition, the Euler-Lagrange equation for H as given by Eq.(p2|) is satisfied, 
then £{h) = P in Eq.(||). 

Let us examine Eq.(p8|) a little more closely. We note that there are three 
conservation laws, whereas there is only one shape equation. To resolve this 
discrepancy, we decompose the space vector f a into its tangential and normal parts, 

f a = r b e h + / a n. (29) 

Note that, in general, the surface tensor f ab need not be symmetric in its indices. The 
surface covariant divergence of f a gives 

V a f a - (V a f " + K\f a ) e b + (V Q / a - K ab f ab ) n , (30) 

where we have made use of the Gauss- Weingarten equations for the surface S, Eqs. 
TT|) , ( p2[ ) The surface projections of Eq.(|28|) are therefore given by 

V Q / Q - K ab f ab = £{h) , (31) 

V a f ab + K b a f a = 0. (32) 
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The first equation expresses the Euler-Lagrange derivative £(h) in terms of f ab and 
f a . Using £ = P, the shape equation takes the remarkably simple universal form 

V a f a -K ab f ab = P. (33) 

Note that only the symmetric part of f ab contributes. 

The conservation of the stress tensor along tangents to the membrane (|3^) is 
independent of £ . They provide consistency conditions on the components of the 
stress tensor, telling us how the tangential stress must respond to a given normal 
stress whether or not the shape equation holds. One can think of these equations as 
the Bianchi identities associated with surface reparametrizations. It is also clear that 
they hold separately for each term contributing to H. 

Our treatment so far has been entirely general: the only properties of the 
hamiltonian we have used are its reparametrization and translational invariance. We 
will now evaluate the Euler-Lagrange derivative and the stress tensor for specific 
models. 



4. Soap bubbles 

In order to illustrate the ideas that have been developed in the previous sections, we 
consider first the case of a surface dominated by surface tension, e.g. a soap bubble. 
The hamiltonian is simply proportional to the area 

H = /i f dA, (34) 
is 

where the constant n is the surface tension. This is the simplest hamiltonian one can 
write down for a surface; it depends only on the intrinsic geometry of the surface. The 
normal deformation of this hamiltonian is given by 

5 ± H = n [ dAK®, (35) 

where we have used the familiar expression relating the Lie derivative along the 
normals of the area element of £ to its mean extrinsic curvature K = 'y ab K a i 1 , as 
follows from Eq. (JlOj) , 

5 ± dA = KdA. (36) 

The tangential deformation of this hamiltonian is simply given by Eq. (0), with 
h = (j,. 

In this geometrical language, the shape equation is given by 

fiK = P. (37) 

The surfaces that extremize the hamiltonian have constant mean extrinsic curvature. 
This is a second-order hyperbolic partial differential equation for the embedding 
functions, X(£). To bring the Euler-Lagrange equations into a more conventional 
form, using the Gauss- Weingarten equations we have that K = — n - AX, where 
A is the surface laplacian. The tangential projections of AX vanish identically. We 
can now peel this expression and its tangential counterpart to recover, 

/iAX = -Pn. (38) 

In this model there is no surface term arising from the normal variation, so that 
the quantity introduced in Eq. (|l^) vanishes identically, 5 a [$] = 0, and only the 
tangential variation contributes to the Noether charge. This feature is unique to this 
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hamiltonian. The invariance of the area hamiltonian under euclidean motions gives 
the stress 

f Q = -/xe°, (39) 
which is not only tangential but also isotropic, f ab = — /i7 ab . 

5. The Helfrich hamiltonian 

Let us now consider the problem of extremizing the bending energy of a surface 
given by Eq.(Q) subject to constraints on V, A (and possibly M). To implement 
the constraints on the area and the integrated mean curvature, we construct the 
constrained hamiltonian 

H = aH {2) + f3(M - M ) + fi(A - A ) . (40) 

This is a sum of terms of the form 



H [n) = J dAK n . (41) 

One could consider a hamiltonian with a more general dependence on the K a b, or on 
its derivatives. However, at quadratic order, a term proportional to K ab K ab is locally 
equivalent to K 2 . This is because the Gauss-Codazzi equation @ relates the difference 
K 2 — K ab K ab to the surface scalar curvature 1Z, and the hamiltonian constructed from 
7Z is a topological invariant, as follows from the Gauss-Bonnet theorem. 

As always, the tangential variation of the hamiltonian is straightforward, see Eq. 
([It]). We note that the normal variation of is given by 

5 ± H {n) = [ dA^+^ + nK 71 - 1 {5 ± K)} , (42) 
is 

where we have used Eq. (|3^) in the first term. Using Eqs. (|l0|) and (|l^) for the 
normal variation of the induced metric and the extrinsic curvature, respectively, we 
find 

5±K = -A$ - K ab K ab § . (43) 

Inserting this expression to Eq. ( fi"2| ) , and performing two integrations by parts to collect 
the derivatives of $ in a divergence, we obtain 

S ± H (n) = [ dA {-nAK™- 1 + K^iK 2 - nK ab K ab )} $ 

; [ dA\7 a {K n - 1 \7 a ^>-\7 a K n - 1 <^} . (44) 

We immediately identify the Euler-Lagrange derivative as 

8 (H {n) ) = -nAK 71 - 1 + K n - L (K 2 - nK ab K ab ) . (45) 

Generically, the Euler-Lagrange equation £ = is of second order in derivatives of 
K abl so it is of fourth order in derivatives of the embedding functions X. In the 
exceptional case n = 1 , note that £ — 1Z follows from Eq. (||) which is of second order 
in derivatives of X. 

The shape equation for the model described by Eq.(^0|) takes the form 

- 2aAK + aK(K 2 - 2K ab K ab ) + (31Z + [iK = P . (46) 
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If the membrane possesses a boundary, appropriate boundary conditions in the 
variational principle are identified by examining the divergence appearing in Eq. (Q) . 
There are that $ = 0, which kills both the second term appearing in the divergence 
as well as the contribution from the derivative of $ along C to the term proportional 
to V Q $. The normal derivative l a V a $ remains. Thus we must set l°V a $ = on C. 

We identify the operator £"[<&] introduced in Eq.(|l8|) which corresponds to 
as the 'wronskian': 

Sf n) [$] = -niK 11 - 1 V a $ - ^ a K n - r ) . (47) 

In particular, if $ corresponds to a background translation, we have 

S a (n) = a • S a (n) [n] = -na • [K n ~ l K ab e b - V^^n] , (48) 

where we exploit the Gauss- Weingarten equation (|l^) to simplify the first term. We 
thus have from Eq. (|26|), the general expression for contribution to the stress tensor 
coming from Hi n \ , 

f ( a n) = {nK 11 - 1 K ab - K n j ab )e b - n\7 a K n ^n . (49) 

Unlike the soap hamiltonian proportional to the area, the stress tensor f? % does 
possess a component normal to the surface. If n = 1, however, note that f?u = 
— the corresponding stress is tangential. If V a K = at any point, the normal stress 
vanishes there. Note that the antisymmmetric part of f ab vanishes identically for the 
geometrical hamiltonians we consider. 

The hamiltonian H( 2 ) is invariant under scale transformations. It is easy to see 
that this is equivalent to 

/( a 2 fc )7a b = 0. (50) 

For the complete Helfrich hamiltonian (ff0|), the stress takes the form 

f a = [aK(2K ab - Kj ab ) + !3{K ab - Kj ab ) - /i 7 ah ] e 6 - 2aV a Kn . (51) 

This is the principal result of this paper. In general, the tangential stress f ab is 
neither homogeneous nor isotropic. Indeed, it may vanish in places. We note, in 
particular, that whereas the parameter /i is the thermodynamic tension, determining 
the response of the energy to a change in the membrane area, it is not the mechanical 
surface tension. 



6. Rotations and torque 

Let us now consider an infinitesimal rotation, <5X = b x X. We have $ = b-Xxn and 
$<! = b ■ X x e a , and the variation ( p0| ) of the hamiltonian associated with a rotation 
reduces to 

5H So = b • f dA [8(h) X x n - V a m Q ] , (52) 

where the surface vector m a is given by 

m a = -S a [X x n] - hX x e a . (53) 

The quantity m a is identified as the torque about the origin acting on E. Invariance 
of H^ a under rotations then implies, using £ = P, 

V a m Q = P X x n . (54) 
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We isolate the contribution due to the couple of f a about the origin 

m a = X x f a + s a , (55) 

where 

s a = X x S Q [n] - S a [X x n] . (56) 
Neither s a nor the couple due to P alone are conserved. An immediate consequence 



of Eqs. (|2J), (|5J), ($§ is the relation 

V a s a = f a x e a . (57) 

We emphasize that this relation between s a and f a does not depend on the shape 
equation (]23|); it holds for each term contributing to H. Note that s a does not involve 
derivatives in K ab other than those already contained in f°. 

We can also expand s a analoguosly to Eq.(|2|) with tangential and normal 
projections, s a = s ab eb + s a n. We have 

W a s a - K ab s ab = ^ e ba f ab , (58) 

W a s ab + K b a s a = V7 e QC f a l cb , (59) 

where we have used e a b — V« e a e l) nQ /V7! an< ^ e ^va is the Levi-Civita density. The 
antisymmetric part of f ab , if present, would contribute to Eq.([58|). For models that 
depend on polynomials of the extrinsic curvature, the symmetric part of s ab vanishes, 
and there is no torque about the normal, s a = 0. 

In particular, for the soap bubble the torque is simply that due to the couple of 
f a . The intrinsic torque associated with -ff( n ), however, is generally non- vanishing, 

s^^nif-Vxn. (60) 

Thus for the Helfrich hamiltonian s a is given by 

s° = (2ttiC + ^)e a xn. (61) 



7. Global conservation 

Applying the divergence theorem to Eq.(|28|) with £ — P provides the global statement, 
F(C) = f ds l a i a = P [ dAn. (62) 

JC JT.Q 

Recall that l a denotes the unit normal on C pointing out of Eq. The total force 
exerted on the area element Eo by the enclosed pressure P is balanced by the internal 
forces exerted on the boundary curves C. If P = 0, F(C) vanishes. For spherical 
topology, or for a contractible loop C on a higher genus surface this is also true for 
the individual closed C. On a non-contractiblc loop, however, we have instead that 
F(C) is a non- vanishing constant vector. For a surface with the topology of a torus, 
there are two distinct constant vectors corresponding to the two topologically distinct 
non-contractible circuits. Note that the modulus of this vector can be identified with 
the Casimir of the euclidean group associated with translations. For a surface of genus 
g, there are g + 1 distinct values. 

The corresponding global statement for the total torque acting on any closed loop 

C 

M(C) := [ ds l a m a = P [ dAXxn (63) 
Jc Js 
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follows from Eq.(j54|). 

The total force F on a loop C on a soap film is given by 

F(C) =-n J dsl, (64) 

where 1 = l a e a is the unit normal pointing out of the surface at a given point on C 
treated as a space vector. The corresponding torque acting on a loop C due to the 
enclosed membrane, M(C) is given by 



M(C) = -/i J dsXxl. 



(65) 



The conservation law, Eq.(|28|) has non-trivial consequences which imply global 
constraints on the tangential projections f ab . Let us dot Eq. (^8|) with X. We have 

X • V a f a = V a (X • f a ) - e a • f a = -P X • n . (66) 

We note that e a • ff, = f^. We integrate over £ to obtain the global constraint on the 
trace of f ab , 



dA f a a = -3PV. 



dAX-n 



We have used the representation 

for the enclosed volume. In addition to Eq. (|C 
dAK ab f ab = PA, 



there is the global constraint 



(67) 



(68) 



(69) 



obtained directly by integrating the projection ( |32| ) over the surface. 
For a soap bubble, the integr ability condition Eq.(|67|) gives 

2/iA = 3PV . (70) 

This also follows as a consequence of stationarity of extremal configurations with 
respect to scaling. We have 

HA[\X] - PV[\X] = X 2 /j,A[X] - A 3 P V[X] , (71) 

which is extremized, with A = 1, when Eq.([70|) is satisfied. 

For the Helfrich hamiltonian, the integrability condition Eq.(67) gives the well 
known scaling identity (see |?J Eq.(3.10)) 

I3M + = 3PV . (72) 



Note that the scale invariant H( 2 ) does not contribute. 
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8. Perturbations 

As an application of the conservation law ( |2q ) let us examine briefly perturbations 
about equilibrium. We will evaluate S£ using our knowledge of the stress tensor and 
compare our result with the one obtained using a 'direct' approach j9j. In particular 
the tangential perturbation is simply 6h£ (h) = Q a d a £(h). For the normal variation, 
it is convenient to consider the linearization of 



5±_ V a f a = (S±£(h)) n + £{h) S± n . (73) 

To express this equation in a more useful form we can use i5j_n = — (V a $)e a , together 
with 

S x V Q f a = W a S ± i a - [V a (K*)] f Q , (74) 

where we have used for the normal deformation of the connection, 5±T b a i > = V a (K<&). 
Now, taking projections of Eq. (f73|), we obtain the "linearized Bianchi identity" 

[V a <5±f a ] ■ e 6 - f b a V a (K$) - -(Va$)£(/i) , (75) 

and 

F7 a 6 ± f a ]-n-f a V a (K$) = 6 ± £(h). (76) 

The task of computing the linearization of the Euler-Lagrange derivative is reduced 
by one order: one needs only to compute the linearization of the stress tensor. 
Let us introduce the quantities 

p a = X • e„ , q = X-n. (77) 

Note that p a = V a X 2 /2 is a gradient. The Gauss- Weingarten equations can be cast 
in the alternative form 

V Q p b = - K ab q + ^ab, (78) 
V a q =K abP b . (79) 

Consider an infinitesimal dilatation under which <5X = AX. Then q represents the 
normal component. As such we would expect 5£[q] — inhom, where S£ is the 
linearized Euler-lagrange derivative and by 'inhom' we mean terms which correspond 
to infinitesimal changes in the parameters of the model. We note that under an 
infinitesimal rotation (three parameters), e w e ab h x e a p b satisfies S£ — 0. Together 
with the (three parameter) infinitesimal translation e = a • n they exhaust the zero 
modes of S£ = 0. A membrane with translation symmetry along the z— axis defines 
a unique loop on the orthogonal x — y plane. The infinitesimal deformation of this 
loop induced by a rotation about this axis has e ~ p, where p = X • t and t is 
the tangent vector to loop. Thus, in this case, both projections of X satisfy the 
linearized equation. The consequences in the planar reduction of the Helfrich problem 
are discussed in detail in [p2] . 



9. Symmetries 

The approach we have developed is particularly appropriate when the membrane 
geometry possesses some degree of spatial symmetry. In particular, let us consider 
a membrane that is axially symmetric. Then we can exploit the integrated form of 
the conservation law given by Eq. ( |62[ ) , and balancing forces on a circle of constant 
latitude we obtain a relation between the appropriate components of the stress tensor 



Stresses in lipid membranes 



13 



and the external pressure. When we specialize our considerations to the Helfrich 
hamiltonian ( fiCf ) , this relation reduces to the well known first integral of the Helfrich 
shape equation for axially symmetric configurations (see e.g. |l5|, |2^, [lq]). 

We choose the standard axially symmetric chart on R 3 : p,z,ip. As coordinates 
on the surface we choose the arc-length along the meridians, I, as well as <p. The 
embedding of this geometry in R 3 can then be expressed as 

p = R(l),z = Z(l). (80) 

The line element induced on the surface is 

ds 2 = dl 2 + R(l) 2 dtp 2 . (81) 

The relation 

R' 2 + Z' 2 = 1 (82) 

defines /, where the prime denotes a derivative with respect to arc- length /. A basis 
of tangent vectors for the surface adapted to this coordinate system is then given by 

ei = (R',Z',0) b v = (0,0,1), (83) 

and the outward pointing normal is 

n= (-Z',R',0). (84) 

Let us consider a configuration with spherical topology. We choose the loop C 
to coincide with a circle of fixed latitude, z = constant. The vectors l a = (1,0), and 
e a = (0, R^ 1 ) denote the unit normal and tangent vectors respectively to this circle 
on the surface. We note that 1 = e a l a = e;. 

We now examine the integrated form of the conservation law given by Eq. ( |62| ) . 
We decompose the components of the stress tensor, f ab and f a , with respect to the 
surface basis {l a , e a } as 

r 6 = /±±H 6 + /||||(7 o6 -n 6 ), (85) 

f a =f±l a + f\\e a , (86) 

where all four coefficients are independent of the angle (p. The potential off-diagonal 
term in f ah is zero by axial symmetry. The term that appears on the l.h.s. of Eq. 
©is 

lj a = f±±l + f±n, (87) 
so that, integrated around the loop, and using ds = R(l)d(p, we have 

ds lj a = 27rJ2(/j.j. <I> + f± (n) ) , (88) 

where (•) denotes an average over ip. All but the z-components vanish. We have 

(1) = (0,Z',0), (n) = (0,R',0). (89) 
For the r.h.s. of Eq.(|62|), this gives 

2vrP J dlR(n) = ttR 2 P (0, 1, 0) . (90) 

Therefore Eq. (|62|) reduces to the relation 



/j_L Z' + f ± R' = ^f. (91) 
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We emphasize that this expression does not depend on the details of the model, but 
only on the assumption of axial symmetry and the choice of the particular loop. 
We define now the angle by sin = Z', cos O = R'. Then 
dZ 

tan0=— (92) 
dR 



When we decompose the extrinsic curvature as in Eq. (p5[), we find 

K±± = &, *:|||| = ^, (93) 
and for the mean extrinsic curvature, 

* = 0' + ^. (94) 
For a soap bubble, the relevant components of the stress tensor are, from Eq. 

9b, 



/±± = -M, /± = 0, (95) 

so that the relation (|9l]) reduces to 

sin = -R/R , (96) 

where Rq = 2/i/P. The unique solution which is regular at the poles is a circle of 
radius Rq. 

Let us now consider a membrane described by the Helfrich hamiltonian (f4(i|). As 
follows from Eq. (^) for f a ( n ) we have 

fti=K { 2K^-K) = & + ^)&-*£), 

sin s 



/(2) =-2K' = -2( 0' 



R 

f^=K ±± -K = -sm&/R, 

Therefore using these expressions, together with the ones in (|95|), the relation ( |9l| ) 
gives 

^ I sin©\' sin©. sin©, . ^ 

-2a cos 6 6' + —— +a(6' + — — )(6' — )sin© 

\ R J R R 

-/S— M sm0 = — . (97) 

This coincides with the first integral for the axysymmetric shape equation obtained in 
p5[ , and which is studied e.g. in Q. Our approach should be useful in providing an 
interpretation for an additional constant of integration C that appears in these works. 

Had we taken another loop, for example, the closed loop running along the 
meridan (p = 0) we find that we obtain a not so useful integral identity. 
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10. Remarks 



In this paper, we have exploited a combination of variational principles, and 
conservation laws to examine the physics of two-dimensional surfaces which are 
described, at a mesoscopic level, by an effective hamiltonian which depends locally on 
the surface geometry. In the case of the Helfrich hamiltonian, using this approach, we 
have derived both the stress tensor and the torque, associated with translational and 
rotational invariance, respectively. The shape equation, which determines equilibrium 
configurations, is now identified as one element of a conservation law. 

The shift in focus onto the stresses in the membrane provides both an intuitive 
description of equilibria and a useful calculational tool. We demonstrated this in 
the context of axially symmetric configurations. Using the integrated form of the 
conservation law for the stress tensor, the first integral of the shape equation associated 
with this symmetry was obtained in an surprisingly economical way. 

We have explored various applications of this framework: a membrane with an 
exposed free boundary in Q ; and the adhesion of a membrane onto a substrate in |26[ . 
While both problems have been addressed previously in the axisymmctric context, 
the boundary conditions associated with the geometry makes them very awkward to 
handle when this symmetry is relaxed. We show how our framework is well suited 
not only to identifying these boundary conditions in this case, but also to providing 
a physical interpretation for them. In a forthcoming paper, we will examine in detail 
how to approach perturbation theory in the same framework. 
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